XXX 



A Kelvin-wave cascade on a vortex in superfluid 4 He at a very low temperature 
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A study by computer simulation is reported of the behaviour of a quantized vortex line at a very 
low temperature when there is continuous excitation of low-frequency Kelvin waves. There is no 
dissipation except by phonon radiation at a very high frequency. It is shown that non-linear coupling 
leads to a net flow of energy to higher wavenumbers and to the development of a simple spectrum 
of Kelvin waves that is insensitive to the strength and frequency of the exciting drive. The results 
are likely to be relevant to the decay of turbulence in superfluid 4 He at very low temperatures. 



PACS numbers: 67.40.Vs, 47.32. Cc, 47.37.+q 

It is well known that quantized vortices can be formed 
in superfluid 4 He. Such vortices can support a transverse 
and circularly polarized wave motion (a Kelvin wave), 
with the approximate dispersion relation for a rectilinear 
vortex [l| 
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where k is the quantum of circulation (h/m^), a is the 
vortex core parameter, and c ~ 1. The existence of these 
waves in an inviscid fluid was first discussed as a theoreti- 
cal possibility in the nineteenth century Q , but an exper- 
imental study in a fluid without viscosity had to await the 
discovery of quantized vortices in superfluid 4 He. Kelvin 
waves in uniformly rotating superfluid 4 Hc were first ob- 
served experimentally by Hall and a number of in- 
teresting experimental and theoretical studies have been 
published subsequently; see, for example, Glaberson et al 
[5| on an instability in the presence of an axial flow of the 
normal fluid; and the study of non-linear effects, leading 
to soliton behaviour [|| and to an associated sideband 
instability Q. Non-linear effects remain interesting, and 
an aspect of them that is important in our understand- 
ing of Kelvin waves at very low temperatures, and which 
may be of rather general interest in non- linear dynamics, 
is discussed in this paper. 

At temperatures where there is a significant fraction of 
normal fluid Kelvin waves in superfluid 4 He are damped 
by mutual friction, which is the frictional force exerted 
on a vortex when it moves relative to the normal fluid. 
This paper is concerned with the expected behaviour of 
Kelvin waves at very low temperatures, when damping 
due to mutual friction can be neglected. Under these con- 
ditions Kelvin waves can be damped only by radiation of 
phonons, but the damping is expected to be extremely 
small unless the frequency is very large, typically of 
order 4 GHz (k ~ 2 nm _1 ). Kelvin waves of lower fre- 
quency are essentially undamped. In these circumstances 
Kelvin waves of a particular low frequency can lose en- 



ergy only by non-linear coupling to waves of a different 
frequency. We are led therefore to consider the follow- 
ing situation. Suppose that we have a rectilinear vortex 
of finite length, and that we continuously drive one low- 
frequency Kelvin mode on this vortex. The mode will 
grow in amplitude, until non-linear effects give rise to a 
transfer of energy to other modes, particularly at higher 
frequencies. This process will presumably continue un- 
til modes are excited that have a frequency sufficiently 
high for effective phonon radiation. The aim of this pa- 
per is to predict the details of this process, by means 
largely of computer Simula tions. We shall consider, for 
example, whether there is a steady state, in which energy 
is injected at a low frequency and dissipated by phonon 
emission at a high frequency, and whether there is some 
well-defined and simple spectrum of the Kelvin waves 
existing in this regime. We shall regard this regime as a 
cascade, although we are not sure whether it is a strict 
cascade in the sense that energy is transferred in steps 
to higher wavenumbers. We shall find that a well-defined 
spectrum does seem to exist, that it is simple in form, 
and that it is remarkably insensitive to the amplitude 
and frequency of the drive. 

This behaviour is interesting in its own right. However, 
we ourselves were led to investigate it by an interest in 
the decay of turbulence in superfluid He at very low 
temperatures @, where mutual friction has a negligible 
damping effect on vortex motion. As in classical turbu- 
lence, energy in superfluid turbulence must probably flow 
from larger to smaller length scales, and it has been sug- 
gested that on the smallest scales the relevant motion is 
a Kelvin wave on a vortex with wavenumber greater than 
the inverse vortex spacing. It is then of interest to under- 
stand how energy can flow in a system of Kelvin waves 
towards higher wavenumbers, where it can ultimately be 
dissipated at the highest wavenumber by phonon radi- 
ation. In the context of superfluid turbulence the de- 
tails are likely to be quite complicated (Kelvin waves are 
likely to be generated by vortex reconnections), and it is 
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not yet generally accepted that a Kelvin-wave cascade is 
strictly necessary. Nevertheless it seems likely that such 
a cascade has a significant role. These matters have been 
discussed in a recent review Q , and we plan further dis- 
cussion in forthcoming papers. 

We believe that the behaviour we find may exist also 
in other types of system and may therefore have other 
applications. It should be explained that similar prob- 
lems to that discussed here have been addressed by other 
authors H [HI El HI > but they relate to situations that 
are significantly different; for example, where there is no 
steady state or where the length of vortex line is con- 
strained not to increase. Thus Araki and Tsubota et al 
carried out numerical simulations on an initial con- 
figuration in which a vortex ring approaches a rectilinear 
vortex, Kelvin waves being generated by reconnections 
when collision takes place. There was no steady input of 
energy and no obvious dissipation. The numerical sim- 
ulations of Kivotides et al jlfl| related to an initial con- 
figuration of four vortex rings, Kelvin waves again being 
generated by reconnections when the rings collide. Again 
there was no steady input of energy and no obvious dis- 
sipation. In both these simulations the authors found 
some evidence for the development of a spectrum simi- 
lar to Eq. (6) below, although the situation is obviously 
different from that considered in this paper. The com- 
putational work of Nemirovskii et al 12] did deal with a 
steady state in which waves on a vortex ring were gener- 
ated by low-frequency noise, although the total length of 
line was kept artificially constant; the results suggested 
the existence of a spectrum similar to Eq. (3) below, 
although it did not exhibit the insensitivity to the driv- 
ing conditions that we ourselves find. Generally, then, 
it seems that this other work has not led to such simple 
and clear-cut results as we report in this paper, and, as 
we shall discuss in later papers, it may be less relevant 
to the decay of superfluid turbulence, especially when, 
as is often likely to be the case jjjj, energy flows from a 
large reservoir associated with large-scale quasi-classical 
turbulent motion into motion on a scale less than the 
vortex-line spacing. 

We consider a model system in which the helium is 
contained in the space between two parallel sheets, sep- 
arated by distance £b = 1 cm, with a single, initially 
rectilinear, vortex is stretched between opposite points 
on the two sheets. Kelvin waves can be excited on this 
vortex, and periodic boundary conditions are applied at 
each end. Thus the allowed wavenumbers of the Kelvin 
waves are given by 



where n is an integer (> 0). We imagine that one of 
these modes, with n equal to a small integer no, is con- 
tinuously driven, so that its amplitude tends continuously 
to increase. As the amplitude increases non-linear cou- 
pling to other modes sets in, and we can expect energy to 
flow from the mode uq to other modes, with both larger 



and smaller wavenumbers. Now we introduce a suitably 
strong damping for all modes with n exceeding a large 
critical value n c . This is intended to mimic the effect of 
phonon emission, although, because of inevitable com- 
putational limitations, it is occurring at a much smaller 
frequency. Then we ask whether there is a steady state, 
described by an energy spectrum Ep., in which the en- 
ergy input to the mode no is balanced by dissipation in 
the modes with n > n c . We find that such a steady 
state does seem to exist, and we determine the character 
of the corresponding energy spectrum. We observe no 
reconnections. 

An important feature of the model lies in the fact that 
Kelvin modes with wavenumbers less than 2tt/£b cannot 
be excited, so that energy cannot flow to smaller and 
smaller wavenumbers. The relevance of this feature to 
the decay of superfluid turbulence will be discussed in a 
later paper. 

The simulations are based on the vortex filament 
model, and they are similar to those described by 
Schwarz [l3| and used in more recent work by one of 
the authors 0]. The undisplaced vortex lies along the 
z-eods. Its calculated motion is based on the full Biot- 
Savart law and therefore takes account of both local and 
non-local contributions. The force that drives one mode 
is of the form V pnsm(koz — u>ot), where ko = 2Tmo/iB, 
p is the density of the helium, and luq is related to ko 
by the dispersion relation (1). Damping at the highest 
wavenumber allowed by the resolution of the simulations 
(1/60 cm) is applied by a periodic smoothing process, 
the details of which will be described in a later publi- 
cation; this allows an effective dissipation at the highest 
wavenumber that can adapt to the flux of energy through 
/c-space arising from the drive. 

Fig. 1 shows how the total length of line evolves in 
time after application of the driving force. We see that it 
reaches a steady average value, suggesting the existence 
of a steady state. 

We express our more detailed results in terms of the 
root mean square amplitudes Ck(t) = (CCfc) 1 ^ 2 of the 
Fourier components of the displacement of the vortex. 
Figure 2 shows how these amplitudes develop in time af- 
ter the application of a drive with V = 2.5 x 10~ 5 cm 
s _1 and ko — lOir cm -1 . We see that initially only the 
mode that resonates with the drive is excited. However, 
as time passes, non-linear interactions lead to the excita- 
tion of all other modes. Eventually the spectrum reaches 
a steady state, shown by the blue line, and there is then 
no further change. In this steady state energy is injected 
at a certain rate at the wavenumber ko, and it is dissi- 
pated at the same rate at thehighest wavenumber. For 
large values of k, where the modes form practically a con- 
tinuum, the steady state is observed to have, to a good 
approximation, a spectrum of the simple form 

H = Ai- B x k~\ (3) 

where the dimensionless parameter A is of order unity. 
Figs. 3 and 4 show the effects respectively of increasing 
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FIG. 1: The development in time of the total length of vortex 
line. 
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FIG. 2: Time development of Cfe(^)- The red, green, black, 
orange and blue lines refer respectively to averages over 0-800 
sec, 10000-10800 sec, 20000-20800 sec, 40000-40800 sec and 
140000-140800 sec. The blue line relates to the steady state. 

the drive amplitude V by a factor of ten and of changing 
the drive wavenumber Icq. We see that there is no effect 
on the steady state, within the error of the simulations, at 
least at the higher wavenumbers. The steady state takes 
longer to be established at the lower drive amplitude, 
which suggests that even with a small drive amplitude the 
same steady state would be established after a sufficiently 
large time; but limitations on the time available for a 
simulation have not allowed us to check this suggestion. 

Important questions relate to the relationship between 
the drive amplitude and the power input to the system 
of Kelvin waves. For small times after the drive is first 
established non-linearities in the system are relatively 
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FIG. 3: Steady-state values of Cfe(^) f° r two different drive 
amplitudes. The blue and red lines are for, respectively, V = 
2.5 x 10~ 5 cm sec -1 and V = 2.5 x 10 -4 cm sec -1 . The mauve 
line has the form of Eq. (3). 
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FIG. 4: Steady-state values of (k(t) for drives at three differ- 
ent wavenumbers. The red, green and blue lines refer respec- 
tively to ko — 2-7T cm -1 , fco = 4tt cm -1 and fco = 10-7T cm -1 . 
Again the mauve line has the form of Eq. (3). 



unimportant, and the power input can be calculated from 
the product of the drive amplitude and the amplitude of 
the velocity response at wavenumber fco, taking proper 
account of phase differences. For later times, however, 
the system exhibits very strong non-linear behaviour, 
and this simple technique is no longer applicable. En- 
ergy seems to be injected through processes involving two 
wavevectors, the difference between which is equal to fco- 
Therefore we have not yet been able to calculate how the 
power input varies with the drive amplitude, although 
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we hope to do so later. Careful analysis of the operation 
of damping in our model for large wavenumbers might 
allow us to obtain the rate of dissipation of energy, but 
we have not yet completed this analysis. For the present 
we can only make the reasonable assumption that an in- 
crease in the drive amplitude does increase the power 
input. We can conclude with some confidence therefore 
that the spectrum is insensitive to the drive amplitude, 
the drive frequency, and the power input from the drive. 

The mean energy per unit length of vortex in a mode 
k is related to by the equation 



E k = e K k 2 ( 2 k 



(4) 



where 6k is an effective energy per unit length of vortex, 
given by 
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It follows from Eqs. (3) and (4) that 
E k = Ae K {kt B r X . 



(5) 



(6) 



We conclude that in our model system a steady state 
cascade does develop, that this state is characterized by 
the energy spectrum (6), and that, remarkably, this spec- 
trum is insensitive to the frequency and amplitude of the 
drive and to the power input at the drive frequency. 

It is interesting to ask what physics underlies this re- 
sult. We suggest that, in the steady state and for waves 
of wavenumber of order k, there is a saturation in the 
local amplitude of the Kelvin waves at a value of roughly 
/c -1 ; this arises from a sudden onset of strong non-linear 
effects when the amplitude is of order the wavelength. 
We make the reasonable assumption that the total mean 
square amplitude of the displacement of the vortex is in- 
dependent of the length (£ B ) of the vortex. The spectrum 
(3) must then be proportional to . Our assumption 
about the sudden onset of non-linear effects means that 
the only other parameter on which this spectrum can de- 
pend is k. The form (3) then follows from a dimensional 
argument. We propose now to investigate whether this 
type of behaviour can be found in other forms of wave 



propagation, at least if they have the same types of dis- 
persive and non-linear characteristics (which are known 
to lead to soliton behaviour |(|). 

As far as we are aware, there are as yet no experimen- 
tal results that relate directly to the behaviour of Kelvin 
waves in superfluid 4 He at very low temperatures. Such 
experiments, which might involve the excitation of Kelvin 
waves on the regular array of vortices existing in the uni- 
formly rotating liquid, as in the early experiments of Hall 
H, would be of great interest. Measurements might be 
made of the rate at which energy is transmitted to the 
array from a suitable transducer, and of the rate at which 
energy eventually appears as heat after the Kelvin-wave 
cascade has been established. 

In summary we have reported the results of computer 
simulations of the behaviour of Kelvin waves on a recti- 
linear quantized vortex of finite length in superfluid 4 He 
at a temperature so low that the waves suffer no attenua- 
tion from mutual friction with the normal fluid, the only 
attenuation arising from phonon radiation at a very high 
frequency. The waves are excited by continuousy driving 
the system at a small wavenumber. The amplitude of the 
driven mode increases until non-linear coupling leads to 
a transfer of energy to all other modes. A steady state is 
established, described by a simple energy spectrum, the 
form of which is remarkably insensitive to the strength 
and other details of the drive. 



Acknowledgments 

We are grateful for stimulating conversations and cor- 
respondence on the general background to this work with 
Carlo Barenghi, Demos Kivotides, Boris Svistunov and 
Sergey Nemirovskii. MT is grateful to the Japan Soci- 
ety for the Promotion of Science for financial support 
through the Japan-UK Scientific Cooperative program 
(Joint Research project) and through a Grant-in- Aid for 
Scientific Research (Grant no. 15340122). WFV is grate- 
ful for financial support through a Royal Society Joint 
Project on Experimental and Theoretical Study of Su- 
perfluid Turbulence. 



[1] R. J. Donnelly, Quantized Vortices in Helium II (Cam- 
bridge University Press, Cambridge, 1991). 

[2] W. Thomson, Phil. Mag. 10, 155 (1880). 

[3] H. E. Hall, Proc. Roy. Soc. A245, 546 (1958). 

[4] W. F. Vinen, Phys. Rev B64, 134520 (2001). 

[5] W. I. Glaberson, W. W. Johnson and R. M. Ostermeier, 
Phys. Rev. Lett. 33, 1197 (1974). 

[6] H. Hasimoto, J. Fluid mech. 51, 477 (1972). 

[7] D. C. Samuels and R. J. Donnelly, Phys. Rev. Lett. 64, 
1385 (1990). 

[8] W. F. Vinen and J. J. Niemela, J. Low Temp. Physics, 
128, 167 (2002). 



[9] B. V. Svistunov, Phys. Rev. B52, 3646 (1995). 
[10] D. Kivotides et al, Phys. Rev Lett. 86, 3080 (2001). 
[11] T. Araki and M. Tsubota, J. Low Temp. Physics, 121, 
405 (2000). 

[12] S. K. Nemirovskii, J. Pakleza and W. Poppe, Russian J. 

Eng. Thermophysics, 3, 369 (1993). 
[13] K. W. Schwarz, Phys. Rev B31, 5782 (1985); Phys. Rev. 

B38, 2398 (1988). 
[14] M. Tsubota, T. Araki and S. K. Nemirovskii, Phys. Rev. 

B62, 11751 (2000). 



